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ABSTRACT 

Context. New and existing large-scale radio surveys of the Milky Way at centimetre wavelengths can play an important role in 
uncovering the hundreds of expected but missing supernova remnants in the Galaxy's interstellar medium. We report on the discovery 
of two supernova remnants (SNRs) designated G152.4-2.1 and G190.9-2.2, using Canadian Galactic Plane Survey data. 
Aims. The aims of this paper are, first, to present evidence that favours the classification of both sources as SNRs, and, second, 
to describe basic parameters (integrated flux density, spectrum, and polarization) as well as properties (morphology, line-of-sight 
velocity, distance and physical size) to facilitate and motivate future observations. 

Methods. Spectral and polarization parameters are derived from multiwavelength data from existing radio surveys carried out at 
wavelengths between 6 and 92 cm. In particular for the source Gl 52.4-2.1 we also use new observations at 11 cm done with the 
Eff'elsberg 100 m telescope. The interstellar medium around the discovered sources is analyzed using 1-arcminute line data from 
neutral hydrogen (Hi) and 45-arcsecond ^^CO(J=1^0). 

Results. Gl 52.4-2.1 is a barrel shaped SNR with two opposed radio-bright polarized flanks on the North and South. The remnant, 
which is elongated along the Galactic plane is evolving in a more-or-less uniform medium. G190.9-2.2 is also a shell-type remnant 
with East and West halves elongated perpendicular to the plane, and is evolving within a low-density region bounded by dense neutral 
hydrogen in the North and South, and molecular (^^CO) clouds in the East and West. The integrated radio continuum spectral indices 
are -0.65+0.05 and -0.66+0.05 for G152.4-2.1 and G190.9-2.2 respectively Both SNRs are approximately 1 kpc distant, with 
G152.4-2.1 being larger (32x30 pc in diameter) than G190.9-2.2 (18x16 pc). These two remnants are the lowest surface brightness 
SNRs yet catalogued at I^igHz ^^"^^ ^ m'^ Hz'^ sr'^ 
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1. Introduction 

Supernovae are the predominant source of new energy and ele- 
memtal input into the interstellar medium (ISM). Shock waves 
from supernova blasts carve out hot ionized tunnels in the ISM, 
sweep up and ionize the gas in the neutral ISM, twist and contort 
the Galactic magnetic field, compress nearby interstellar clouds 
and could also trigger star formation and create the soup of cos- 
mic rays in the Galactic environment that plays an important 
role in pressure balance in the Milky Way ISM and star forma- 
tion (among other things). Besides the detailed understanding 
of individual physical properties of supernova remnants (SNRs) 
that is needed to get a more global picture of their place in the 
ISM, a simple accurate count of the number of Galactic SNRs 
is required. Assuming a mean a ge of radio shell-type SNRs of 
>60,000 yrs ( [Frail et al.|[T994l ) and a rate of one SNe per 30 - 
50 years in spiral galaxies like the Milky Way (IKoo et al. 2006)), 
the ISM should play host to between one and two-thousand radio 



supernova remnants at any given epoch. However in 2009 only 
274 SNRs were catalogued in the Milky Way ( [Green||2Q09l ), and 
although this has increased to >31Q today in light of rece nt ob- 
servations (see cataloguj^of [Ferrand & Safi-Harb 2012) there 
are still many "missing" SNRs, undoubtedly due to difficulties 
in identifying very faint objects. The "missing SNRs problem" 
is also apparent in models of the angular distribution of H ii re- 
gions and SNRs in the Galaxy which predict about 1000 Galactic 
SNRs dLietal.l [19911 ). 

Many radio-emitting SNRs descend from type lb, Ic and 
II (core collapse) supernova events, and are more likely to be 
evolving in a low density environment created by the stellar 
wind(s) of their massive and powerful progenitor (and cluster 
companions). The majority of Galactic SNRs in the catalogue 
of |Kothes et al. ( 2QQ6 |) are found in the adiabatic expansion 
phase, where energy from their shocks is being efficiently con- 
verted into radio. X-ray and/or optical radiation. This indicates 
an observational bias favouring SNRs that are significantly in- 
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teracting with their local ISM environment, and suggests that 
many SNRs, whose shocks are still more-or-less freely expand- 
ing inside pre-made cavities and have not encountered the in- 
ner walls of the bubbles created by stellar winds, remain unob- 
served. Discoveries of new SNRs within high- sensitivity, high- 
resolution and wide-field radio continuum and polarimetric sur- 
veys of the Galactic plane (for example Brogan et al. 2QQ6[ 
Gao et aL) 2011} are becoming common now that such sur- 
veys are more widely available. This success at alleviating the 
missing SNRs problem is being used to motivate a new gen- 
eration of surveys (e.g. GALFACTS, POSSUM, LOFAR, and 
surveys to be performed by the SKA). Their sensitivity to all an- 
gular structures and their wide fields allows detection and study 
of large, low-surface brightness objects in the context of their 
local Galactic background, while their high-resolution lowers 
the confusion limit, allowing one to distinguish small angular 
structures from large ones (and to separate them). The Canadian 
Galactic Plane Survey (CGPS, iTaylc^ et al.|[2Q03i ) is the original 
such survey, and probes the ionized, neutral, and magneto-ionic 
ISM over unprecedented spatial dynamic range. Since 2001 this 

dataset has enabled the disc overy of 7 new SNRs in the 1 ^^ and 
■)nd 



2^^ quadrants of longitude (Kothes et al. 2001 



Kothes 



2003 



[Kothes et al.) |20Q5l [Tian"et al.| |2007| [Kerton et al.| POOV), and 
allowed new insights into the SNR population in our Galaxy. In 
this paper we add two of the faintest-known SNRs to the list of 
catalogued remnants in our Galaxy. 

The new SNR candidates in this paper were discovered in 
two individual CGPS 21 cm mosaics that had been source- 
subtracted and subsequently smoothed to gain S/N and better re- 
veal large unseen extended emission in the survey. The two low 
surface-brightness shell-type SNRs are designated G152.4-2.1 
and G190.9-2.2 (names from Galactic coordinates of their 
geometric centres). Flux densities at five radio wavelengths 
(/16 cm, 11, 21, 74 and 92 cm; in frequency v =4.8 GHz, 
2.7 GHz, 1420 MHz, 408 MHz and 327 MHz respectively) are 
measured with Stokes I data. Polarized intensity maps with B- 
field vectors from Stokes Q and U maps are presented at 
11, and 21 cm (for G190.9-2.2 only 6 cm and 21 cm polar- 
ization data are available). Steep integrated radio spectra for 
each of the respective 1?5 and 1° diameter shells of Gl 52.4-2.1 
and G190.9-2.2 show the predominantly non-thermal nature of 
the sources. Finally we present 1-arcminute resolution Hi and 
^^CO(J=1^0) line channel maps towards each remnant and esti- 
mate LSR velocities, kinematic distances and physical diameters 
for each. 

The primary purpose of this paper is to announce the dis- 
covery of the objects, present evidence for their identification 
as SNRs and provide only the basic measurements that can be 
reliably gleaned from survey data. It is hoped this will aid fu- 
ture specialized observations of both objects, observations with 
deeper sensitivities which will be needed to embark on detailed 
studies of their physics (e.g. spatial and frequency variation of 
spectral indexes, relative amount of free-free thermal and sy- 
chrotron emission, age and shock velocities, etc), beyond the 
scope of this discovery paper. 

2. Observations 

2.1. 21 cm and 74 cm 

The 21 cm Stokes I, Q, U continuum and H i line, and 74 cm con- 
tinuum data come from the CGPS, which is described in detail 
by [Taylor et aL| ( |2003] ) The release of th e CGPS 21 cm Q and U 
maps are described by Landecker et al. (2010j). For our purpose 



(a qualitative look at the polarization structure of the objects) we 
use the interferometer-only Q and U data without the large-scale 
structures restored. CGPS data were observed with the 7-antenna 
Synthesis Telescope (ST) of the Dominion Radio Astrophysical 
Observatory (DRAG). Near G152.4-2.1 the elliptical synthe- 
sized beam of this interferometer in the continuum is 66V 3x49^.^6 
at 21 cm, and 3:8x2:8 at 74 cm. Near G190.9-2.2 these values 
become highly elliptical at 2:6xO:83 and 9: 1x2:9. The practice 
of observing all baselines from 617mtol2.9m over 144 hours 
per field leads to good sensitivity despite the relatively small 9-m 
diameter of the primary antennae: in the vicinity of G152.4-2.1 
(G 190.9-2.2) measured 1-sigma sensitivities are ATg =60 mK 
(40 mK) in the 21 cm continuum, and 0.8 K (0.4 K) at 74 cm. 
In 21 cm Stokes Q they are -100 mK (42 mK) and in Stokes 
U are 80 mK (41 mK). The more elongated synthesized beam 
near £ =190° accounts for its lower noise. Emission from struc- 
tures larger than those the shortest baseline of the ST can sample 
are added in from surveys by the Eff'elsberg 100 m, the Stockert 
25 m and the DRAG 26 m telescopes. Data from single-dishes 
like these are of particular importance to restoring the true flux 
density of faint extended objects like new SNRs, which is missed 
by the incomplete spatial frequency sampling of ST images. The 
CGPS data used here are two of 84 mosaics designated "MSTl" 
and "MEQl", which are all available to the community through 
the Canadian Astronomy Data Centre (CADCQ). 

2.2. 1 1 cm 

11 cm data of the two new SNRs were obtained with the 
Efl'elsberg 100-m telescope. For G152.4-2.1 new observations 
were made between August 2006 and January 2007 with a two 
channel cooled receiver with a bandwidth of 80 MHz centered 
at 2639 MHz. The two circularly polarized 80 MHz wide chan- 
nels were correlated by an Intermediate Frequency polarimeter, 
which records Stokes I, Q and U simultaneously. In addition, 
the 80 MHz bandwidth was split into 8x10 MHz wide chan- 
nels to eliminate narrow-band interference and to measure possi- 
ble large rotation measure (RM) variations. 3C286 served as the 
main calibrator with a flux density of 10.4 Jy, 9.9% of linear po- 
larization at a polarization angle of 33°. We observed a 2°Ax2°A 
large field centred on G152.4-2.1 by moving the telescope beam 
either along Galactic longitude or along Galactic latitude direc- 
tion with a speed of 37sec. The individual scans of each map 
have T separation providing full sampling of the 4:4 beam of the 
Eff'elsberg telescope at 2639 MHz. The standard data reduction 
software package, based on the N0D2 format ( Haslam 19741 
was used for the Eff'elsberg continuum observations. The maps 
in NOD2-format were edited for interference, and "scanning ef- 
fects" were suppressed by applying the unsharp masking method 
( [Sofue & Reich |1979| ). We ended up with nearly 15 coverages 
which were combined into the final map by using the weaving 
method described by Emerson & Graeve ( 1988 ). The rms-noise 
of the final maps was measured to be 2.5 mK (~ 1 mJy) for Stokes 
I and 1.5 mK for Stokes Q and U. Analysis of the narrow chan- 
nel polarization data does not indicate the existence of large RM 
features exceeding 100 rad m"^ or more. Earlier low resolution 
RM data by S poelstra (198 4 ) indicate RMs around rad/m^ in 
this direction of sky. By defining the pixels at the end of each 
scan to be zero, large-scale emission exceeding the size of the 
map may be lost, which is not an issue for studying a discrete 
continuum source within a field, but might aff'ect Stokes Q and 
U. 
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For G190.9-2.2 we extracted total intensity data from the 
EfFelsberg 11 cm (2695 MHz) Galactic plane surve}[^ de- 
scribed in I Furst et ([1990)), which has a sensitivity of about 
ATg =16 mK. For details of the receiver used for this survey 
and the antenna response of the Effelsberg 100-m telescope see 
[Reich et al.| ( [T984l ). 

2.3. 6 cm 

Data at A6 cm (y =4.8 GHz) in Stokes I, Q and U are from 
the Sino-German polarization survey of the Galactic plane, made 
with the Urumqi 25 m telescop e at Nanshan Station, China (Sun 
et al.| 2007 [ Gao et al. 2010). Basic parameters of this survey 
are a resolution of 9^5, and a measured 1-sigma sensitivity in 
our data of ATg ~1.5 mK in I and ~0.5 mK in Q and U. We 
have used the original Urumqi polarization data without the ab- 
solute WMAP corrections. These data already include emission 
of several degrees in extent (much larger than our objects) and 
so are appropriate for polarization studies of all but the largest 
discrete objects. 



2.4. 92 cm and Other Frequencies 

Stokes I data for G152.4-2.1 were taken from the Westerbork 
Northern Sky Surve>0 (WENSS, [Rengelink et al.| |l997t) at 



Table 1. Integrated flux and spectral properties of G152.4-2.1 
and G190.9-2.2 from point-source-subtracted Stokes I maps at 
five radio frequencies. 



92 cm (325 MHz). As this survey only extends north of decli- 
nation -h28.5 degrees, data for G190.9-2.2 is unavailable. The 
measured 1-sigma flux density variation in these data is about 
3.5 mJy/beam, or ATg ^10 K with the synthesized beam of 
54''x54'7sin(^). Due to the restriction on short observing base- 
lines for the Westerbork interferometer (about ~36 m, or 1.5 
times an antenna diameter), the WENSS survey misses flux from 
structures greater than about 1?4 degrees in size. This will impact 
the 92 cm flux density derived for G152.4-2.1 only in a small 
way. 

A faint exte nded so urce at the position of Gl 90.9-2.2 was 
catalogued by Kassim| ([1988 ) at 30.9 MHz. He did not iden- 
tify the object he catalogued as NEK 190.9-2.3 as a SNR. 
Nonetheless the orientation and size of emission contours in his 
maps are similar to G190.9-2.2, and it is likely these are one and 
the same object. The integrated flux density listed by Kassim for 
NEK 190.9-2.3 is S^^^y^^ =13.6 Jy+20%. 

Dust emission from the objects was searched for in IRAS 
12, 25, 60 and 100 yum, and MSX 8 jim infrared survey data, but 
no emission spatially correlated to either source's radio morphol- 
ogy was found at any of these infrared wavelengths. The com- 
bined Wisconsin H- Alpha Mapper (WHAM) and Virgina Tech 
Spectral line Survey (VTSS) map made by |Finkbeiner ( 2003| ) 
was also searched for Ha emission from each object (angular 
resolution -6'). We describe the discovery of emission corre- 
sponding with the radio appearance of G190.9-2.2 in Sec|T2 



3. Appearance, Integrated Flux Densities and 
Spectra 

To better represent extended continuum emission, unresolved 
"point" sources were first removed from each map. For 21, 74 
and 92 cm data, sources down to Signal/Noise- 2(T were readily 
found by automated source finder "findsrc", which examines an 
image using a matched "point-source" wavelet filter to enhance 



Flux (+AS) 


G152.4-2.1 


G190.9-2.2 


S4812[Jy] 


1.19+0.15 


0.49+0.07 


S2639/2695 [Jy] 


2.07+0.28 


0.57+0.10 


Sl420[Jy] 


3.16+0.75 


1.21+0.20 


S408[Jy] 


5.85+1.39 


1.74+0.51 


S327[Jy] 


7.18+1.61 




a (Soc y^) 


-0.65+0.05 


-0.66+0.05 


Sv=lGHz[Jy] , ^ ^ 


3.54+1.81 


1.30+0.63 


2 1 GHz [ 1 0-23 Wm-^Hz- 1 sr" ^ ] 


5.32+2.71 


4.11+2.00 


Centre (^, b) 


152?57,-2?04 


190P97,-2n3 


Size X %ij^xangle 


99:6 X 94:8x0° 


69:3x60:4x30° 



point-like sources and locate them in an image. Sources are mod- 
elled with two-dimensional elliptical Gaussians plus a base level 
using the program "fluxfit". Both routines are part of the DRAG 
Export Software Package, a modular software package of gen- 
eral and telescope-specific programs for processing and analysis 
of images made by the DRAG Synthesis Telescope ( [Higgs et aT 



J997 ). 6 and 1 1 cm maps had relatively fewer unresolved sources 
since their resolution is much poorer, so the brightest individ- 
ual sources were identified using the original 1-arcminute 21 cm 
map overlaid as a guide, and these were modelled and subtracted 
with fluxfit. For brevity we do not list derived source fluxes and 
coordinates in this paper. To aid in identifying the shell and its 
edges the 21, 74 and 92 cm maps were next smoothed to 4:4 (the 
resolution of the 11 cm maps). Noise for these final smoothed 
maps are 25 mK, 0.4 K, and 1.6 K (at 21, 74 and 92 cm respec- 
tively) for G152.4-2.1, while the noise level for G190.9-2.2 is 
15 mK and 0.4 K at 21 and 74 cm respectively. It is these source- 
removed and smoothed maps that are shown in Figures [T] and [2] 
White 21 cm contours are overlaid on each map for comparison, 
also derived from the smoothed 4:4 resolution map. 

3.1. G1 52.4-2.1 

We refer here to total intensity maps at /16 , 11, 21, 74, and 
92 centimetres (4.8 GHz, 2.7 GHz, 1420 MHz, 408 MHz and 
327 MHz) in Figure [T] All total (I) and pol arized intensity (PI) 



maps are presented in the colour-scheme of |Green| ( |2011 ). The 



^ Available at http ://w w w3 . mpif r-bonn . m pg . de/survey.html | 
http://www.astron.nl/wow/testcode.php 



emission structure of G152.4-2.1 is best seen in the 21 cm map. 
The very bright H ii region Sh2-206 is seen in the top-right cor- 
ner of I maps (and 6 cm PI map). The radio emission from this 
SNR candidate consists of a pair of enhanced shells of emission 
with steep outer edges in the North and South. This is the typ- 
ical "double loop" or "barre l"-shaped structure of a pure shell- 
type supernova remnant (see [Manc hester 1987[ |Gaensler||1998[ 
for example), where two shells of enhanced radio emission are 
found atop a more diff'use central emission plateau which is 
likely generated by those parts of the source that are expand- 
ing towards us or moving away from us.'Gaensler (1998 ) found 
that the orientation of these barrel- shaped SNRs is mostly par- 
allel to the Galactic plane - as is the case for G152.4-2.1. He 
suggests that the reason for this shape is a combination of the 
uniform magnetic field structure which is preferentially paral- 
lel to the Galactic plane and the expansion into ISM structures 
which are elongated in longitudinal direction. [Caswell (19771 
also describes asymmetries in shell remnants as the outcome of 
ISM density gradients. 
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We determine the polarization properties of G152.4-2.1 us- 
ing Stokes Q and U polarization maps at 6, 11, and 21 cm. 
In Fig. [T] we display polarized intensity images with B -vectors 
overlaid. The length of each vector is determined by the polar- 
ized intensity at that point; however, as the intention is simply 
to highlight the field's direction, the scale in cm per mK is arbi- 
trary from map to map in Figures [T] and |2] Synchrotron emission 
from supernova remnants should be inherently polarized up to 
70%; however Faraday rotation by the ISM and internal eff'ects 
will depolarize this signal, and this depolarization is much more 
significant at 21 cm than at 6 cm, eff'ectively scrambling low- 
frequency polarized intensity from a SNR. This eff'ect can be 
easily seen in our polarized emission images. While there is a 
strong correlation between total power and polarization at 6 cm, 
the ambient polarization features become more and more domi- 
nant at longer wavelength. 

G152.4-2.1 shows polarized intensity at 6 cm well corre- 
lated to continuum features, a good indicator that the emission is 
synchrotron in nature. The fractional polarization on the bright- 
ened regions of the shell is between 50 and 60%. G152.4-2.1 
also seems to sit on an elongated slightly curved feature in po- 
larized intensity more than 2 degrees in length. This large po- 
larization filament seems to have a tangential intrinsic magnetic 
field as indicated by the B -vectors. At this high frequency, where 
eff'ects of Faraday rotation should be small, the B -vectors of the 
SNR candidate indicate a tangential magnetic field in the north- 
ern shell, which is expected for a mature shell-type SNR. The 
southern shell shows a slight rotation of up to 30 degrees clock- 
wise from a tangential structure, which requires a rotation mea- 
sure of about -130 radm"^. This would produce a rotation of 
about 90 degrees for polarization angles at 1 1 cm and we find the 
B -vectors in the 11 cm polarization image in this area to be al- 
most radial. We did not produce a rotation measure map because 
at the lower frequencies the background/foreground polarization 
structures become more and more dominant and it is difficult to 
distinguish between them and genuine SNR emission. 

From the polarized emission structure we conclude that the 
SNR is likely interacting with the large polarization filament, 
the nature of which is unclear. Its bright polarized emission im- 
plies an enhanced, compressed internal magnetic field. In this 
scenario the progenitor star of Gl 52.4-2.1 would have exploded 
at the upper edge of this filament, sweeping up its magnetic field 
explaining the filament's gap at the location of the SNR at 6 cm 
and the similarity of the B -vectors of both objects. 

The 21 cm PI map is shown in Figure [T] The only morpho- 
logical correlation between the PI and I maps at 21 cm seems to 
be the distinct dark channel of depolarization that snakes West to 
East across most of the South shell's face, which tightly follows 
the continuum around the South-East corner of G152.4-2.1. It is 
moderately consistent with a compressed magnetic field aligned 
in the East- West direction along the shell which cancels with the 
principally North-South field seen in the unrelated bright PI from 
the ISM immediately above it. Some of this bright PI emission 
on the shell's Western face at i =152?2,/? = -2?4 could originate 
with the SNR's shell, but how much is impossible to say from 
these maps, as is its precise role in creating the depolarization 
"snake". 



3.2. G1 90.9-2.2 

Total Power maps at 6, 11, 21 and 74 centimetres (4.8 GHz, 
2.7 GHz, 1420 MHz, and 408 MHz) are shown in Figure [2] 
The first impression of G190.9-2.2 is one of an incomplete 
shell of emission with a brightened and sharpened Southern 



edge to the object. On closer inspection however we believe 
that Gl 90.9-2.2 shows a barrel-shaped structure very similar 
to G152.4-2.1, with two lobes of enhanced shell emission on 
either side (East and West) of more diffxise emission in the cen- 
tre. The Eastern lobe is seen well at 21 cm and 11 cm, and the 
Western half is bright at 6 cm. The angle of the barrel's line of 
cylindrical symmetry on the sky is about 110 degrees from the 
Galactic longitude axis. If there is a correlation between this ori- 
entation of the "barrel" and the ambient magnetic fiel d and e lon- 
gation of ambient structures as proposed by G aensler ( |1998 ), we 
must have either a slight anomaly here where the local magnetic 
field has a large angle with Galactic longitude or G 190.9-2.2 is 
a SNR expanding inside a stellar wind bubble. Since the radio 
surface brightness of this SNR candidate is very low it is likely 
expanding in a low density environment, and in fact G190.9-2.2 
does appear within a gas-poor circular region bounded by both 
dense molecular clouds in the East and West, and dense neutral 
hydrogen caps in the North and South (see Section]?]). 

Optical emission atop G190.9-2.2 is seen in WHAM-hVTSS 
maps: a bright isolated unresolved patch of Ha emission some 
~5 Rayleighs above background is found to correspond with the 
brightest part of the 21 cm emission (the South-East portion of 
the radio shell). The patch of Ha is shown in the 21 cm panel 
of Figure [2j outlined with dotted contours. The origin of this Hor 
emission is clearly seen in an optical (red) POSS-II image in 
Figure |3] the outline of which is shown as the dashed box atop 
the 21 cm panel in Fig.|2] The patch corresponds to a very long, 
thin filament of red nebulosity that runs parallel to the 200 mK 
21 cm contour (the innermost black contour around the brightest 
inner portion of the South-East shell on the 21 cm panel of Fig.|2j 
the same contours are drawn in white on Fig.|3]). This likely is 
highly compressed cooling post-shock gas just downstream of 
the SNR shock wave, which is expanding in a South-Easterly 
direction. 

The PI map at 6 cm (Fig. [2]) shows the eastern shell to be 30- 
50% polarized, but the other parts of the SNR candidate seem 
to be unpolarized. This could be the result of the low surface 
brightness combined with uncertain contributions from the unre- 
lated interstellar medium along the line-of-sight. The B -vectors 
derived for the eastern shell indicate a tangential magnetic field 
which would be expected for a mature shell-type remnant. This 
provides further evidence that Gl 90.9-2.2 is a supernova rem- 
nant. 

3.3. Integrated Flux Densities 

Integrated flux densities are summarized in Table [T] Tabulated 
fluxes are averages obtained over 20 polygons each with two 
kinds of background surfaces fitted to the data values at the 
polygon vertices: twisted-plane and twisted-quadratic surfaces. 
Twisted plane surfaces have the property that orthogonal cuts 
through them at one particular rotation angle are linear, and 
twisted quadratics have one orientation in which all cuts through 
the surface are linear, and all orthogonal cuts are quadratic. 
Uncertainties in fluxes (AS) accumulate from three sources: i) 
the uncertainty in the total intensity calibration of each map 
(±4, 10, 5, 15 and 20% for A =6, 11, 21, 74, and 92 cm respec- 
tively), ii) the deviation in the background- subtracted flux over 
difl'erent polygons enclosing the emission (drawn by eye), and 
iii) the mean difl'erence in background flux levels estimated with 
both types of background surface fits. A linear fit of log S ver- 
sus logy was calculated with least squares weighted by AS. Each 
object shows a steep power-law radio spectrum identifying the 
integrated emission as non- thermal (see Fig. [4]). For G152.4-2.1 
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Fig. 3. Cutout of the bright radio shell of G190.9-2.2 corre- 
sponding to the box outlined on the 21 cm panel of Fig.[2j show- 
ing optical (red) emission from the Second Palomar Observatory 
Sky Survey (bright emission is shaded dark here). A thin com- 
pressed filament of red nebulosity is seen running parallel and 
next to the innermost 21 cm radio contour, which is the com- 
pletely enclosed white contour in the centre (the three white con- 
tours here are the same as are in the 21 cm panel of Figure [2]). 
This filament likely traces shocked gas downstream of the SNR 
shock. 



the integrated spectral index is a = -0.65+0.05 (Sy oc y^). For 
G 190.9-2 .2 a = -0.66+0.05 (without the 30.9 MHz flux point 
of |Kassim| ( p^88l ) a = -0.59+0.14). Also reported in Table [T] 
is each object's estimated flux density and surface brightness at 
1 GHz, the mean centroid coordinates, and the major and minor 
angular diameters and orientation angle of each shell (fitted by 
eye). At 1 GHz the surface-brightness of both SNRs (as listed in 
Table [T]) would make them the faintest yet catalogued, below the 
previous faintest S NR G156.2+5.7 (5.8x10""' ^" ' 
[Reich etal.|p^92l ). 



4. H I and CO Velocities and Distances 

We determine systemic LSR velocities for both SNRs along their 
lines-of-sight using -T resolution CGPS 21 cm Hi line data 
and 45'' resolution 2.6 mm (1 15.3 GHz) i2cO(J=1^0) line data 
from the Exeter FCRAO CO Galactic Plane Survey (described 
in |Mottram & Brunt|[20T0l [Brunt et aL| 2012). 

Figure |5] is a 9-channel montage of ^^CO and Hi towards 
G152.4-2.1. Each channel is spaced to 0.82 km s"^ and the 
velocity resolution is 1.32 km s"^ The CO maps have been 
smoothed to V spatial resolution to match the H i. In all velocity 
channels, line emission towards G152.4-2.1 shows no obvious 
cavity nor shell of neutral or molecular gas surrounding it (as 
an SNR evolving in a stellar wind bubble might; e.g. 3C434.1, 
Foster et aL| [2004[ ). Rather, the continuum boundary is sharply 
delimited by a wall of bright CO that systematically wraps the 
Eastern boundary of the SNR from the South to the North. This 
bright CO (T^ =1-2 K) emission begins outside of the Southern 
shell boundary at vlsr = -8.9 km s"^ and in each of the fol- 
lowing channels (-10.5 to -13 km s"^) the CO neatly matches 
the curvature of the sharp South-East comer of the shell. Also 
tracing the SNR's South-East and East periphery are numerous 



Fig. 4. Integrated radio flux spectra for each object, and overlaid 
power-law fitted with least squares weighted by uncertainties. 
For Gl 90.9-2.2 the spectrum with index a4 is calculated over 
408 MHz-4.8 GHz (A14 -6 cm) only and does not include the 
30.9 MHz flux point from [Kassim[ ( [1988| . 



Hi Self- Absorption (HIS A) clouds (white contours in Fig. [5]) in 
CGPS maps at the same velocities. The CO and HIS A environ- 
ment seems to form a "wall" of cool dense gas on the sky next 
to the East continuum boundary of the SNR. As well, CO and 
HIS A appear to neatly fill in an Eastern "bay" formed by the 
continuum contours at ^ =152?9, = -1?9. An isolated bright 
cloud of H I emission (single 90 K thin black contour in each 
map of figure), the brightest such cloud in the entire field and 
the only H i emission feature possibly associated with the SNR, 
is seen exactly in the gap between the CO clouds that wrap along 
the outside of the Southern shefl, at ^ =152?2, b = -2?6. AH of 
these spatial coincidences occur in nine channels between -8.9 
and -15.5 km s"^ Their spatial relation with respect to the shell 
is similar in channels -10.5 to -13 km s"^ suggesting a sys- 
temic velocity of vlsr - -12+2 km s"^ for G152.4-2.1. 

Figure [6] is a 9-channel montage of ^^CO and Hi cen- 
tred on G190.9-2.2. Two bright CO clouds that begin to ap- 
pear at +6.0 km s"^ flank the SNR on its left (East) and right 
(West) sides. This relationship is particularly striking in chan- 
nels from +5.1 to +3.5 km s"^ where the CO clouds are dense 
and bright and continue to flank the SNR's East and West contin- 
uum boundary. At the same velocities an extensive HISA cloud 
that well traces the outline of the CO cloud in the East ap- 
pears, and another HISA cloud just begins to appear atop the 
Western CO at +3.5 km s"^ In later channels vlsr <2.7 km s"^ 
this CO+HISA appears atop the SNR's face, suggesting it is in 
front or behind it and sporting little relation with the contin- 
uum appearance. Hi emission surrounding G190.9-2.2 is much 
more complex than seen towards G152.4-2.1, and appears in 
the same velocity channels as the CO but in the opposite cor- 
ners: above and below the object's North and South boundaries 
respectively. Both of these North and South Hi "caps" follow 
the curvature of the SNR's limb at 21 cm in four channels +7.6 
- +5.1 km s"^; thereafter the North arc disappears and the South 
has grown into a large H i wall ofl' of the South South- West cor- 
ner of the SNR. The overall appearance of Gl 90.9-2.2 being 
nestled within CO+HISA walls to the East and West and Hi 
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emission shells to its North and South is particularly striking in 
the channel map at vlsr = +5.1+1.6 km s"^ Based on this ap- 
pearance we assign this as the systemic velocity of G190.9-2.2. 

Even with a systemic velocity in hand, ascribing distances to 
these two objects is not straightforward. Under undisturbed cir- 
cular rotation vlsr = -12+2 km s"^ places G152.4-2.1 within 
the "Local arm" at a distance of -1.1+0.1 kpc (using the rotation 



curve of Foster & Cooper , 2010 ). This SNR then has a linear di- 
ameter of 32x30 pc. The gradient Ar/AvisR here is 110 pc per 
km s~^ , which gives the minimum uncertainty. However the asso- 
ciation with HIS A implies that neutral hydrogen at two different 
distances is present in these velocity channels (i.e. a cold fore- 
ground and a warm background), showing that circular rotation 
in this direction is anything but undisturbed and that the error on 
the distance is expected to be greater. 

A kinematic distance to G190.9-2.2 of 1.0+0.3 kpc is indi- 
cated based on the CO+Hi-based velocity of +5.1 km s"\ im- 
plying physical dimensions of 18x16 pc. The uncertainty in this 
distance is large however, as the 1/sin^ projection here causes a 
small velocity range to represent a large line-of-sight path, and 
small velocity deviations from circular rotation are magnified 
into large distance uncertainties (here Ar/AvisR ~200 pc per 
km s"^). We know such deviations are present as remarkable 
and extended HISA is seen through the field, so again the dis- 
tance error is a lower limit. 



5. Conclusion and Future 

The purpose of this paper is twofold. Firstly we introduce two 
extended faint discrete objects discovered in the COPS and show 
evidence (mainly through their radio spectral and polarization 
properties) that classifies them as SNRs. Secondly we provide 
a basic qualitative interpretation and quantitative physical prop- 
erties (fluxes, distances, sizes) that will aid the community in 
planning future new observations of them at radio and espe- 
cially other wavelengths (e.g. X-ray, optical). The objects were 
discovered by deep point source subtraction (down to 2cr level) 
and subsequent smoothing of two selected CGPS 21 cm con- 
tinuum mosaics, to increase the S/N of faint extended emis- 
sion. Their identity as SNRs was then discovered spectrally 
through comparison with other wavelengths (e.g. 74 cm com- 
ponent of the CGPS; the Sino-German 6 cm polarization sur- 
vey). Undoubtedly many more such extended low surface bright- 
ness SNRs await discovery within this rich dataset using this ap- 
proach. Such discoveries will further address the "missing SNR" 
problem that is a key driver of new and planned deep radio sur- 
veys of the Milky Way's ISM. 
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Fig. 1. Total (I) and polarized intensity (PI) contoured maps of Gl 52.4-2.1. The emission features (colour bars in brightness tem- 
perature units of K, at right) are themselves outlined in thin black contours (levels indicated in the text above each panel). B -field 
vectors are overlaid on PI maps (length proportional to PI; scale diff'ers between maps). Resolutions are 9' 5 at 6 cm and 4f4 at 11, 
21, 74 and 92 cm (except 21 cm PI which is V). PI maps are missing large-scale backgrounds. Four 21 cm contours (white) at 
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11cm I + contours (30,60,90mK) 74cm I + contours ( 1 . 1 , 1 .5, 1 .9,2.3K) 
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Fig. 2. Total (I) and polarized intensity (PI) contoured maps of Gl 90.9-2.2. The emission features (colour bars in brightness tem- 
perature units of K, at right) are themselves outlined in thin black contours (levels indicated in the text above each panel). B -field 
vectors are overlaid on PI maps (length proportional to PI; scale diff'ers between maps). Resolutions are 9f5 (6 cm), 4f3 (11 cm, 
21 cm I), 2!6x0'85 (21 cm PI), and 9f9 (74 cm). The box outline on the 21 cm panel corresponds to the optical image in Figure |3} 
the elliptical patch of Ha emission identified in WHAM-hVTSS data is also shown (dotted; 11,12 Rayleighs) on this panel. Three 
21 cm contours (white) at 100,150, and 200 mK are shown overlaid on all maps, except the 21 cm panel which has white contours 
drawn from the 9f9 74 cm I map. 
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Fig. 5. ^^CO line channel maps towards G152.4-2. 1 . The brightness scale is antenna temperature T^ in K, and the angular resolution 
of the maps is T. Three thick white-on-black contours traced at 100, 130 and 160 mK delineate the total power appearance of 
G152.4-2.1 at 21 cm. Overlaid on each map are thick pure white contours that outline HISA clouds seen in CGPS Hi maps at the 
same velocities. A single Tg =90 K contour (very thin black) outlines the brightest Hi cloud in the field, appearing in between the 
molecular clouds that follow the outside of the southern shell. All of this activity peaks slIvlsr = -12 km s"^. 



9 



Foster et al.: Two New Supernova Remnants 




Fig. 6. ^^CO line channel maps towards Gl 90.9-2.2. The brightness scale is antenna temperature T^ in K, and the resolution V. 
Three thick white-on-black contours traced at 100, 150 and 200 mK delineate the total power appearance of G152.4-2.1 at 21 cm. 
Extended HISA is contoured with thick pure white lines (level -20 K below the mean Tg level of each H i channel), and H i emission 
at 6 and 9 Kelvin (above the mean of each map) is shown with the thin black contours. In channels with visr = +4.3 and 5.1 km s"^ 
molecular gas is seen bracketing the SNR's East and West boundaries, and curved shell-like atomic gas clouds bracket the North 
and South boundaries. 
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